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A detailed view of a Molecular Cloud in the far outer disk of M 33 
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ABSTRACT 

The amount of Ht present in spiral galaxies remains uncertain, particularly in the dim outer regions and in low-metallicity envi- 
ronments. We present high-resolution CO(I~0) observations with the Plateau de Bure interferometer of the most distant molecu- 
lar cloud in the local group galaxy M 33. The cloud is a single entity rather than a set of smaller clouds within the broad beam 
of the original single-dish observations. The interferometer and single-dish fluxes are very similar and the line widths are in- 
distinguishable, despite the difference in beamsize. At a spatial resolution of 10 pc, beyond the optical radius of the M 33, the 
CO brightness temperature reaches 2.4 Kelvins. A virial mass estimate for the cloud yields a mass of 4.3 x lO** and a ratio 
A'(H2)//co(i-0) - 3.5x 10^" cm"-/(K km s"'). While no velocity gradient is seen where the emission is strong, the velocity is redshifted 
to the extreme SW and blue-shifted to the far NE. If the orientation of the cloud is along the plane of the disk (i.e. not perpendicular), 
then these velocities correspond to slow infall or accretion. The rather modest infall rate would be about 2 x lO^'^Mgyr"'. 
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1. Introduction 

How much molecular gas is present in spiral galaxies? What de- 
termines whether or not a molecular cloud will form in a given 
place and time? Do all molecular clouds form stars and, if so, at 
what rate? These are central questions to understanding where 
and why stars form in galaxies yet they are difficult to answer. 
One of the major difficulties is the absence of direct observa- 
tion of cool H2. The standard tracer of molecular gas, the CO 
molecule, which is the most abundant molecule with a rota- 
tional dipole moment, is clearly a poor tracer of H2 in very low- 
metallicity environments dRubio et alJll99ll:lBolatto et al.ll2008h 
and contested i n other situations, particularly in the cool outer 
parts of spirals ( iPfenniger et al.lll994tlSodroski et alj|1995h . CO 
emission has been detected in only a very few galaxies be- 



^ ■ yond the nomina l optical radius R25 (e.g. |D igel et al. 1991 



X 



iBrand et al.lll987) for the Galaxy and lBraine & Herpia (,20041) 
and Braine et alJ ( 2007h for other galaxies. 

The Local Group spiral M 33, R25 = 30.8' ~ 7.5 kpc, is 
an ideal candidate to better study the above issues because it is 
very nearby, such that molecular clouds can be detected individ- 
ually, with a classical spiral disk structure with no geometrical 
un certainties, and an on ly slightly subsolar metallicity (-factor 
2, iMagrini et al.l |2010|) . The observation s presented h ere are 
part of a series by our group dGardan eFal. 2()07: Gratier et alJ 
120101 120 121: iKramer et al.l l201 0: Brai ne et al.ll2010bir on the gas 



and du s t of M 33 with s imilar goals to work by e . g. WilsonI 



(Il997h . lEngargiola eTaP (l2003h . iRosolowskv et all (l2007h . or 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



iBigiel etalJ ( l2010h . In this work, we present high resolution ob- 
servations with the Plateau de Bure interferometer of the very 
distant cloud detected as M33-18 bv lBraine et al.l (1201 Oah . just 
beyond R25 and the most distant detected to date. This is roughly 
equivalent to twice the solar circle distance. 

Being able to spatially resolve clouds is important to deter- 
mine whether the weak emission in the outer disks of spirals is 
due to the very low brightness temperatures of the clouds or to 
their scarcity in a telescope beam encompassing an area much 
greater than that of a molecular cloud. Increasing the resolution 
enables us to study the cloud itself to determine whether it is a 
single cloud and to estimate the mass necessary to gravitation- 
ally bind it together. This is the goal of the cuiTent work. 



2. Observations and data reduction 

The most distant molecular cloud to date in the disk of M 33 has 
now been observed at high resolution in the CO(l-O) line. The 
data were acquired in early 2009 with the interferometer on the 
Plateau de Bure operated by the Institut de RadioAstronomie 
Millimetrique in the 'C and 'D' configurations. A single field 
was observed centered on RA (J2000) 01:34:12.5 and Dec 
(J2000) 31:10:32. 

Calibration was straightforward and standard routines within 
CLIC were used. 3C454.3 and MWC 349 were observed once at 
the beginning of each track respectively as bandpass and flux 
flux calibrators, 0234-)-285 was used as phase and amplitude cal- 
ibrator and was observed every 12 minutes for the remainder of 
the track. The calibration uncertainties at 115 GHz at the PdBI 
are estimated to be 10%. 
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Datacubes were made at 4 different spatial resolutions 
by varying the weighting scheme applied to the visibilities. 
Tapering at 50 and 90 meters with natural weighting yielded a 
beam of HPBW 6.65" x 5.46" (PA 99°) and 4.50" x 3.49" (PA 
114°) respectively. Pure natural weighting with no taper, in or- 
der to maximize point source sensitivity, produced a beam of 
HPBW 3.58" X 2.67" (PA 119°) and robust weighting gave a 
resolution of 2.74" x 2.40" (PA 105°). These angular resolutions 
correspond to linear resol utions of 10 - 25 pc for an assumed dis- 
tance to M 33 of 840 kpc dGalleti et alJ2004 . The initial spectral 
resolution was 0.1 km s ' (39 kHz) but in all cases two chan- 
nels were combined such that the final channel separation was 
0.2 km s"'. Each cube has 76 channels. All image processing, 
including deconvolution, was done in the MAPPING environ- 
ment (see http://www.iram.fr/IRAMFR/GILDAS). 

Interferometric observations always have a complex beam 
pattern so deconvolution ("cleaning") is necessary and is often 
the most delicate part of the data reduction. The dataset with 
the lowest angular resolution has the highest brightness sensitiv- 
ity and is thus most sensitive to extended structures. The "dirty" 
cubes produced as per the previous paragraph were first cleaned 
with the default clean box, which covers much of the central re- 
gion and yields an unbiased indication of where signal is present. 
For the 22 channels which appeared to show emission, polygons 
were drawn to indicate the regions to clean and this enabled 
us to further improve the imaging. The Clark clean algorithm 
was used with a maximum of 10000 iterations and 500 major 
cycles except for one channel where we used 1000 iterations. 
The loop gain in all cases was 0.1. The higher resolution cubes 
were cleaned in the same way but the polygons from the low- 
resolution data were used for the cleaning in order to avoid the 
bias of following positive noise at high angular resolution. 



3. Morphology, dynamics, and mass of the distant 
outer disk cloud 

As can be seen in Fig.[Tl the cloud presented here is much fur- 
ther out in the disk of M 33 than others detected up to now. 
Such clouds are rare and difficult to find (see also cloud M33- 
20 in iBraine et al ] l2010ah . Figure |2] shows the integrated inten- 
sity (moment zero) maps at the four resolutions. We have drawn 
what we consider to be the cloud outline on the panels in Fig. [2] 
except for the upper left panel which shows the lowest resolu- 
tion - but highest brightness sensitivity - data where we draw 
a contour including possible extended emission, a sort of maxi- 
mum cloud size derived from the sum of the polygons used for 
cleaning the data. Table [1] provides the integrated fluxes within 
the se contours. 

iBraine et a D (l20T0al) measured an integrated intensity of 
0.61 K km s or 2.9 Jy km s in the CO(l-O) line at 21" reso- 
lution with the IRAM 30 meter telescope but the pointing center 
was 3 arcseconds south of the cloud peak (which was not known 
before these observations). As can be seen from Table 1, very 
little if any flux is resolved out by the present interferometric 
measurements. The 'cloud' contour in Table 1 has an area of half 
that of the 30meter telescope beam and the 'extended' emission 
contour is about 1/3 larger than the 30m beamsize. 

Figure [3] shows the structure of the source as a function of 
velocity, integrating over about 1 km s"' and covering the veloc- 
ity extent of the CO emission. While almost all of the flux is 
contained within a 3 km s"' interval, there is a separate patch of 
weak emission at lower rotation velocity detected to the south- 
west, seen in the lower right panel of Fig. [3] 



Figure |4] shows the spectra at four positions in the cloud 
at the four resolutions. The full line width at zero intensity is 
4.4 km s"' and at half power the line width is 2.4 - 2.6 km s"' 
measured directly on the integrated spectrum. The line is 12 or 
13 channels wide, depending on the exact region over which 
the spectrum is summed. Interestingly, the line width does not 
change with the angular resolution - the highly tapered data have 
the same linewidths as the high resolution cubes. Also, the line 
width changes little with position over the cloud, at least with 
the ~ 10 pc resolution presented here. The line width measured 
with the 30 meter te lescope is the same: Av = 2.5 + 0.2 km s"' 
(IBraine et al.ll20I0al) . 

Even at the highest angular resolution, the cloud does not 
break down into smaller distinct entities, with the possible ex- 
ception of the "low-vel" region seen to the south-west. The 
cloud covers about about 250 square arcsec or 4000 pc^. Unlike 
IBraine et al.l (|2010a), h ere we are able t o measure the size of the 
cloud. Following iSolomon et al.l (Il987l) . we equate R,, - y/A/n, 
and S sm - ■\/7tR^/3.4 - "sfA/dA. Let us note however, that 
while ISolomon et al.l (1 1 9871) delimited clouds at a brightness of 
several Kelvin, the cloud size we measure extends close to the 
zero-intensity level. Adopting their parameters for consistency. 



3fp^/Ao-l/3AG. With f,, = 2.9 and cr^ = Av/2.35, 



we 



M,, 

estimate a mass of 4.26 x 10" * Mp. This is higher b y about 30% 
than the previous estimate bv lBraine et al. l (l2OT0al) based on the 
single-dish spectra. The flux within the 4000 pc^ region in the 
robust map is 3.18 Jy km s"'. The conversion ratio 



-^(H2)//cO(l-0) - 



2m„D^-Q. c^-SAw 



where the first term on the right is the H2 column density and 
the second term is I^^; D is the distance to M 33, nip is the 
proton mass, //, is the hydrogen fraction by mass (0.73), Q. is 
the solid angle occupied by the cloud, and S Av is the velocity 
integrated flux. For the robust map flux of 5 Av =3 . 1 8 Jy km s ' , 
A^(H2)//co(i-o) = 3.5 X 102"cm-V(Kkm s-i). 

If instead we use the rather uncertain extended emission con- 
tour, then the mass becomes 7.07 x 10^ M© and A^(H2)//co(i-0) = 
4.4 X 10^" cm~-/(K k m s~' ). Our estimat es of A^(H2)//co are at 
the upper end of what lLerov et al.l (|20I ll) estimated for M 33. 

Table 1. Fluxes for the extended emission and cloud contours 
for each angular resolution. The first figure in each column is 
the flux measured before primary beam correction and the sec- 
ond after correction. We estimate the statistical uncertainties (in- 
cluding the uncertainty in the cloud contour) to be ~ 15%. 



taperSO 

Jy km s"' 



taper90 
Jy km s"' 



natural 

Jy km s"' 



robust 
Jy km s"' 



cloud 
extended 



2.5, 
3.4, 



2.7 
3.9 



2.8, 3.0 
3.6, 4.1 



2.9, 
3.6, 



3.1 
4.1 



3.0, 
3.7, 



3.2 
4.2 



In the highest resolution map, not tapered and with robust 
weighting (lower right panel of Fig. 2), the peak brightness tem- 
perature reaches 2.4 K. Although not extreme, a brightness tem- 
perature of 2.4 K at 10 pc resolution is actually quite high, par- 
ticularly considering the weak radiation field surrounding the 
cloud. The emission in the GALEX Far-UV band just to the East 
of the cloud (see Fig. |5]l suggests a Far-UV field several times 
the value in the solar neighborhood, only some of which will be 
felt by the cloud, depending on the solid angle occupied by the 
cloud seen from the stars. 
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Figure |5] shows the CO emission in comparison with other 
wavebands. While IR or Ha emission at or near the position 
of the cloud is present, it is typically 2 orders of magnitude 
wea ker than towards star-forming regions in the disk of M 33 
(e.g. lMookerjea etal ]l20Tll iBraine itaPIIoTl although the lat- 
ter have only slight ly stronger CO emission. Using Eq. 9 of 
ICalzetti et al.l (l2007h . the 24 /im, which is the waveband which 
most clearly shows star formation associated with the molecular 
cloud, yields a star formation rate (SFR) of 4 Myr"' . This is 
equivalent to a conversion rate of gas into stars of M(H2)/ SFR 
~ 10"^ years, close to an order of magnitu de slower than in spiral 
galaxies in ge neral (e.g. 'Kennicutt"1998') and M 33 in particular 
jGardan et a Li i,2007: Gratier et al...201Q) . 

The CO(2-l) emission of this cloud, while not m apped, is 
twice as strong in brightness units ( iBraine et al. l l20T0ah towards 
the center; given the cloud size, this implies that the gas kinetic 
temperature and density are sufficient to efficient ly excite the 
J = 2 - 1 transition (see e.g. discussion in App. 1 of lBraine et alj 
Il993h . Since little or no flux is resolved out by the interferome- 
ter, we can smooth to the ~ 11" beam of the 30 meter telescope 
in CO(2-l); the observations presented here confirm a CO(2- 
1)/CO(1-0) ratio close to unity at a scale of 1 1". The neighbor- 
ing UV emission sources cannot heat the interior of the cloud 
due to the high optical depth. The level of star formation, as es- 
timated using the 24 or 8 //m emission, is much lower than in 
many other clouds in M 33 with similar CO emission. 

Ho w does this extreme o uter disk cloud compare with other 
GMC5 .'^ lOratier et al.1 (1201 2l) observed a large sample of GMCs 
with an angular resolution of 12" and found consistently larger 
line widths. This is not due to the angul ar resolution , however, 
because the line width at 21 " resolution ([Braine et al . 2010a) is 
the same as that found here. iGratier et al. ( 12012 ) compared the 
M 33 cloud sample with Milky Way clouds and found that the 
M 33 clouds had smaller line widths than Milky Way inner disk 
clouds - such that the line width of our outer di sk cloud (M33- 
18) is far below that of the lSolomon et alj (Il987h sample of first 
quadrant molecular clouds, in a comp letely empty region in their 
Fig. 1. Perhaps more surprisingly, the lBigiel et alj ( l2010h clouds 
in M 33 (considerably closer to the center, near RA 01:34, Dec 
30:55) and Large Magellanic Cloud sample a lso have broader 
lines, whether for star-form ing clouds or not dKawamura et alj 
120091) . On the other hand, flie lDigel et an(ll994 sample of Milky 
Way outer disk clouds is very similar, with slightly narrower 
lin es and somewhat s maller sizes. The very d istant cloud 2 found 
bv lDigel etaP (11994 and studied further bv lRuffle etaT](l2007h 
is similar in CO luminosity and line width to M33-18 although 
smaller both in Virial mass and siz e, using the average of the dis- 
tances used bv lRuffle et all (l2007h . [Brand & WouterlootI (Il995h 
also find that the Milky Way outer disk clouds have linewidths 
about a factor 2 smaller than inner disk clouds of the same size. 
The CO brightness at lOp c resolution is also similar to far out er 
Galaxy molecular clouds dBrand et alJll987HDigel et al.l[T994l) . 

4. Gas inflow onto the molecular cloud 

Although the H i emission at the position of the cloud is not par- 
ticularly strong, the cloud is part of an arm-like structure seen in 
Hi and it follows the same orientation. Looking at Fig. [3] more 
closely, we see that the blue-shifted emission is to the NE and 
the red-shifted region to the SW. Assuming that the spiral arms 
are trailing, the West side of M 33 is the near side and thus, as 
long as the cloud extensions are in the plane of the galaxy (rea- 
sonable given that they follow the arm orientation), the blue and 
red-shifted emission both correspond to infall towards the cen- 



ter of the molecular cloud. Figure [6| shows the position- velocity 
(PV) diagram for the cloud, following the general orientation of 
the HI arm (see Fig. 5 in Gratier et al. 2010), which is roughly 
NE-SW. The infall velocity is of order 1 km s which, although 
supersonic and thus capable of heating the gas, is not expected to 
create shocks strong enough to destroy grain mantles. The large- 
scale velocity gradient due to the rotation curve is negligible. 

Large-scale infall h as thus far only been detected in DR21 
( [Schneider et al.l 1201 Ol) although i nfall towards a numb er of 
cloud cores has been detected (e.g. lChavarria et al]|2010l) . typ- 
ically via inverse P-Cygni profiles. If we really are tracing infall 
at the scale of 10 - 20 pc, we can estimate the accretion rate as a 
check. By measuring the emission to the SW and NE in respec- 
tively the lower right and upper left panels of Fig. [3] we find that 
fluxes of 0.2 Jy km s"' and 0. 15 Jy km s"' come from the appar- 
ently infalling gas, representing about 10% of the total CO(l-O) 
flux from the cloud. Assuming that the CO emission allows us 
to trace H2 mass, some 4300 M© may be falling onto the central 
GMC. For a distance of 20 pc and a velocity of 1 km s the 
resulting rate is 4300 Mo/2 x lO^yr ^ 0.0002 Mg/yr. Such an 
accretion rate is highly plausible and an order of magnitude less 
than what was estimated for DR21. The presence of two inde- 
pendent clouds cannot be totally excluded but they must either 
be much further apart than their projected distance or must have 
highly differing velocities perpendicular to the line of sight (i.e. 
not observable). 

Averaged over the estimated cloud area of about 4000 pc", 
the cloud surface density is really very low, some 11 Mopc"^. 
However, this is an order of magnitude greater than the aver- 
age stellar surface density so far from the center of M 33, so 
the presence of the GMC represents a real local center of grav- 
ity. Presumably, this is what enables such weak inflow to be de- 
tected. 
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Fig. 1. North-east ern part of M 33 in the IRAC 8 /im band 
(IVerlev et al.ll2007l) . The cloud we present here is in the small 
circle at the top center of the image, representing the primary 
beam of the Plateau de Bure interferometer at the frequency 
of CO Q-O). The "Lonely Cloud", a fair ly CO-strong inter-arm 
GMC dGardanll200"7tlBraine et"al1l2010al) . is also indicated with 
a small circle. The small rectangle at Dec. 30:52 and the small 
triangles forming an 'X' in and around the rectangle indicate the 
individual positions of the [C ii] HIFI observations of the H ii re- 
gion BCLMP 691 presented in Braine et al (2012). The large 
lower rectangl e shows the region aroun d BCLMP 302 observed 
with PACS by iMookeriea etal] (1201 Ih . The center of M 33 is 
at the position of the triangle to the lower right of the image of 
M 33. The giant Hii region NGC 604 can be seen as the very 
bright region to the left around Dec. 30:47 and is marked with a 

X. 
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Fig. 2. Morphology of the outer disk cloud; shown are the mo- 
ment zero or integrated intensity maps made with the tapered, 
natural, and robust weighting. The beam sizes are indicated as 
ellipses in the lower left corner of each map. The color scales 
are shown to the right for the tapered and untapered images sep- 
arately and are in units of Jy/beam. The extended emission con- 
tour is shown in the upper left box and the cloud contour is drawn 
in the other panels. These two contours are used to compute the 
fluxes given in Table [T] 




Fig. 3. Velocity structure of the outer disk cloud at the resolution 
of the map tapered at 90 meters. The cross marks RA 01 :34: 12.3, 
Dec 31:10:35 and is what we later refer to as the (0,0) position; 
the triangles in the panels indicate the high-velocity, peak, and 
low-velocity components shown in Fig.|4l 
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-870 -265 -260 -270 -265 



Velocity (km/s) 

Fig. 4. Spectra as a function of position in the cloud. The (0,0) 
position refers to the cross in Fig. |3] The "peak" position cor- 
responds to the triangle in the upper right panel of Fig. [3] The 
"high-vel" position corresponds to the triangle in the upper left 
panel of Fig. [3] The "low-vel" position corresponds to the trian- 
gle in the lower right panel of Fig.|3] The y-scale varies because 
the spectra are expressed in Jy/beam so the larger beams have 
greater fluxes. The tapered and untapered spectra are shown in 
the same order as in Fig.|2] 




Fig. 6. PV-diagram of the outer disk cloud along a cut with a po- 
sition angle of 49° East of North, made with the datacube tapered 
at 90meters. The Y-axis indicates recession velocity in km s"' 
from slower rotation (-258 km s"') to faster rotation velocities. 
The X-axis indicates the distance from the center in arcseconds, 
with negative values being to the south-west. The color scale is 
in Jy per beam. The bar to the lower left indicates 4", as the beam 
measures 4.50 x 3.49". 
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Fig.5. Comparison of the CO(l-O) emission at the resolution of the natural weighting - contours at 0.05, 0.1, 0.2, 0.3 Jy km s ' - 
with the emission in the Far-IR, Ha (arbitrary units), Far-UV (arbitrary units), H i, and SyumPAH bands. From left to right (top row) 
units are: MJy/sr, Jy/pixel, arbitrary, MJy/sr and (bottom) Jy/pixel, MJy/sr, arbitrary, and integrated Hi line intensity in K km s"'. 
Labels indicating the wavebands are provided in the upper left corners of each pan el and the color wedge is shown above. The H i 
data is at the full 5" resolution of the VLA mosaic described in lGratier et al.l (l2010l) . 



